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mation of two-dimensional unsteady airfoil theory. Six coupled mode shapes, calculated from free 
vibration about the equilibrium operating condition, are used in the linearized stability analysis. 
The study emphasizes the effects of two types of structural coupling that strongly Influence the 
stability of hingeless rotor blades. The first structural coupling is the linear coupling between 
flap and lead-lag bending of the rotor blade. A structural coupling parameter is Introduced 

to simulate variations in flap-lag structural coupling that arise for blades having nonuniform 
stiffness distributions. The second structural coupling is a nonlinear coupling between flap 
bending, lead-lag bending, and torsion deflections. Results are obtained for a wide variety of 
hingeless rotor configurations and operating conditions in order to provide a reasonably complete 
picture of hingeless rotor blade stability characteristics. The stability of torsionally flexible 
blades is strongly influenced by the effects of the bending-torsion sturctural coupling. Without 
precone, typical configurations are usually stable except for low values of (R or low torsion 
frequencies. Addition of precone is strongly destabilizing for a wide range of configurations. 
Except for very low torsion frequencies, the results also Indicate that the structural terms in 
the torsion equation dominate the torsion inertia and damping terms which permits the use of an 
approximate, but simplified, system of equations with fewer degrees of freedom. Finally, the accu- 
racy of the results is sensitive to the number and type of mode shapes used in the analysis. 
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Ames Research Center 
and 

Ames Directorate 

U.S. Army Air Mobility R&D Laboratory 


SUMMARY 


The stability of elastic flap bending, lead-lag bending, and torsion of 
uniform, untwisted, cantilever rotor blades without chordwise offsets between 
the elastic, mass, tension, and aerodynamic center axes is investigated for 
the hovering flight condition. The equations of motion are obtained by simpli- 
f lying the general, nonlinear, partial differential equations of motion of an 
elastic rotating cantilever blade. The equations are adapted for a linearized 
stability analysis in the hovering flight condition by prescribing aerodynamic 
forces, applying Galerkin’s method, and linearizing the resulting ordinary 
differential equations about the equilibrium operating condition. The aero- 
dynamic forces are obtained from strip theory based on a quasi-steady approxi- 
mation of two-dimensional unsteady airfoil theory. Six coupled mode shapes, 
calculated from free vibration about the equilibrium operating condition, are 
used in the linearized stability analysis. The study emphasizes the effects 
of two types of structural coupling that strongly influence the stability of 
hingeless rotor blades. The first structural coupling is the linear coupling 
between flap and lead-lag bending of the rotor blade. A structural coupling 
parameter (H is introduced to simulate variations in flap-lag structural 
coupling that arise for blades having nonuniform stiffness distributions. The 
second structural coupling is a nonlinear coupling between flap bending, lead- 
lag bending, and torsion deflections. Results are obtained for a wide variety 
of hingeless rotor configurations and operating conditions in order to provide 
a reasonably complete picture of hingeless rotor blade stability characteris- 
tics. The stability of torsionally flexible blades is strongly influenced by 
the effects of the bending-torsion structural coupling. Without precone, 
typical configurations are usually stable except for low values of (H or low 
torsion frequencies. Addition of precone is strongly destabilizing for a wide 
range of configurations. Except for very low torsion frequencies, the results 
also indicate that the structural terms in the torsion equation dominate the 
torsion inertia and damping terms which permits the use of an approximate, 
but simplified, system of equations with fewer degrees of freedom. Finally, 
the accuracy of the results is sensitive to the number and type of mode shapes 
used in the analysis. 


*Ames Directorate, U.S. Army Air Mobility R&D Laboratory 



INTRODUCTION 


The general problem of helicopter aeroelastic stability involves coupling 
between the motion of the individual blades and coupling between the rotor and 
the body of the helicopter; control system dynamics may also be involved. The 
complexity of the general problem poses a considerable challenge to the ana- 
lyst, both in developing an analytical model of the system and in understand- 
ing its physical behavior. An important part of the general rotor-body dynamic 
system is the single blade rotating about an axis fixed in space. For many 
problems of practical interest, blade-to-blade and rotor-body couplings are 
not significant and the analysis of a single rotor blade constitutes an impor- 
tant problem by itself. Even when coupling with other blades and the body is 
significant, the single blade behavior usually remains recognizable and can be 
helpful in understanding the behavior of the more complete system. For this 
reason, the dynamics of a single blade forms an important fundamental building 
block in the study of helicopter dynamics. 

Helicopter rotors with cantilever blades are commonly termed hingeless 
rotors. In contrast with the more conventional articulated rotor, the canti- 
lever blades of the hingeless rotor are attached directly to the hub without 
flap or lead-lag hinges. This configuration reduces mechanical complexity and 
improves helicopter flying qualities by increasing rotor control power and 
angular rate damping. The lack of hinge articulation also alters the struc- 
tural characteristics of the rotor blade and can significantly influence aero- 
elastic stability. 

Aeroelastic stability of the articulated rotor blade is primarily deter- 
mined by the relative chordwise position of the aerodynamic center, center of 
mass, and elastic axes. When these are unfavorable, an instability may occur 
that is usually characterized by coupled flap bending and torsion deflections 
with a frequency near the torsion natural frequency. 

In the case of the hingeless rotor blade, another type of instability is 
possible because of the structural coupling between bending and torsion deflec- 
tions of cantilever blades. This type of instability is usually characterized 
by coupled flap bending, lead-lag bending, and torsion deflections of the 
blade with a frequency near the lead-lag bending natural frequency. The struc- 
tural coupling of cantilever blades is significantly dependent on the specific 
configuration parameters of the rotor blade. The magnitude and variability of 
this coupling make the analysis of cantilever rotor blades a complex and impor- 
tant subject. Moreover, the stability characteristics of hingeless rotor 
blades are still not well understood, even though considerable work has been 
devoted to this problem recently- The present investigation is intended to 
provide an improved understanding of the basic structural aspects of these 
problems; therefore, the effects of chordwise offsets between the elastic, 
mass, tension, and aerodynamic center axes are not included in this report. 

The general equations of motion for a torsionally flexible cantilever 
rotor blade contain important nonlinear inertial and structural terms. These 
nonlinearities may be expected to cause nonlinear blade motion behavior, such 
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as limit cycles or nonlinear instabilities. For the present report, this type 
of behavior is not considered, and only the linear stability characteristics 
are treated. This does not mean, however, that nonlinear terms in the equa- 
tions are ignored. In the present analysis, perturbation equations are 
obtained by linearizing the nonlinear equations about a suitable equilibrium 
condition. The influence of the nonlinear terms is then manifest in the 
coefficients of the linear perturbation equations. 

If it is necessary to treat the complete nonlinear equations, more elab- 
orate analyses must be used such as numerical integration of the equations or 
asymptotic expansion techniques. However, these methods are more difficult to 
apply than the standard eigenvalue methods used for linear equations and they 
are not as well suited to investigations where extensive numerical results are 
desired for a broad range of system parameters. Furthermore, an analysis of 
the nonlinearities cannot be justified logically until the linear behavior of 
the system is sufficiently well understood. For these reasons the scope of 
the present report is restricted to a linear stability analysis based on 
linearization of the general nonlinear equations of motion found in reference 1. 
Before briefly describing the procedure used in this report, it will be useful 
to review previous related research. 


Discussion of Previous Research 

Recent investigations of f lap-lag-torsion stability of elastic canti- 
lever rotor blades were preceded by simpler analyses of flap-lag stability 
of torsionally rigid blades. These flap-lag investigations revealed some 
of the complexities of nonlinear aerodynamic and inertial flap-lag coupling 
terms of cantilever rotor blades. The early flap-lag analyses usually 
relied on an approximate representation consisting of a rigid blade with 
spring restrained hinges at the hub to simulate bending flexibility. At 
that time, the most comprehensive equations applicable to rotating elastic 
beams were the linear equations developed by Houbolt and Brooks in refer- 
ence 2. 

In the early flap- lag analysis of reference 3, Young investigated the 
possibility that nonlinear aerodynamic and inertial coupling terms could be 
destabilizing for hingeless rotors. Subsequently, Hohenemser and Heaton 
(ref. 4) showed that linearized equations of motion for the spring restrained, 
hinged, rigid blade approximation could be used to investigate flap-lag insta- 
bility in hover and forward flight. Further investigations in hover by 
Ormiston and Hodges (ref. 5), again using the rigid blade approximation, 
showed that the distribution of flexibility inboard and outboard of the pitch 
bearing (the degree of flap-lag structural, or elastic, coupling) was an 
important factor influencing stability. Flap-lag stability of elastic canti- 
lever blades with uniform properties was studied by Hodges and Ormiston 
(refs. 5 and 6), based on a derivation of nonlinear partial differential equa- 
tions suitable for elastic hingeless rotor blades. The results confirmed that 
the approximate rigid blade model could accurately predict lead-lag damping 
provided that the flap-lag structural coupling was properly represented. 
Similar equations were studied by Friedmann and Tong (refs. 7 and 8), except 
that flap-lag structural coupling was not included. The stability 
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characteristics were determined by using an asymptotic expansion procedure 
with multiple time scales. The results indicated that certain configurations 
shown to be unstable from a simple linearized stability analysis would exhibit 
large amplitude limit-cycle behavior with the more elaborate analysis. 

Concurrently with the restricted flap-lag stability analyses, efforts 
were also made to investigate the complete problem including torsional blade 
deflections. Torsion deflections of hingeless rotor blades are strongly influ- 
enced by nonlinear structural moments caused by flap and lead-lag bending. 

This bending-torsion structural coupling, discussed by Mil’ et at, (ref. 9) is 
proportional to the product of the flap and lead-lag bending curvatures and 
the difference between the two bending flexibilities. Other investigators 
have also noted the importance of bending-torsion structural coupling 
(refs- 10-13) and it has been approximated in some stability analyses using 
the rigid, hinged blade approximation (refs. 14-16). The latter analyses also 
include fuselage degrees of freedom as well as rotor blade deflections. 

One of the first efforts toward the nonlinear flap-lag-torsion problem 
for a fully elastic blade was a derivation of nonlinear equations of motion 
for elastic flap bending, lead-lag bending, and torsion of rotating beams by 
Arcidiacono (ref. 17). The equations were derived in a reasonably detailed 
and complete manner; however, some small nonlinear inertial terms were not 
included - 

Friedmann and Tong (ref. 7) and Friedmann (ref. 18) also developed flap- 
lag-torsion equations. Elastic flap and lead-lag bending of a uniform blade 
were considered, but the torsional deflections were approximated by rigid body 
pitching motion (root torsion) . The results of reference 7 indicated that 
torsion motion was important and that the stability characteristics were sen- 
sitive to the number and type of assumed bending mode shapes used. Flap-lag 
structural coupling was not included in these reports, however, and only a few 
numerical results were presented. 

Hodges (ref. 19) developed a system of nonlinear flap-lag-torsion equa- 
tions for an elastic torsionally flexible cantilever blade and applied them to 
a cantilever blade stability analysis. The equations contained some small 
erroneous structural terms later shown to be due to incorrect shear strain- 
displacement relations. Numerical results illustrated the importance of elas- 
tic torsional flexibility and the possibility of omitting, for some configura- 
tions, torsional dynamics from the equations of motion and retaining only the 
structural (kinematic) effects of torsion. 

Hodges and Dowell (ref. 1) corrected and extended Hodges’ derivation of 
the nonlinear flap-lag-torsion equations using both Hamilton’s principle and 
the Newtonian method. Prior to formal publication of reference 1, Hodges and 
Ormiston (ref. 20) applied these equations to a uniform rotor blade of solid 
cross section. Without blade precone, most hingeless rotor configurations 
with uniform blade properties and typical bending and torsion frequencies were 
found to be stable. Positive precone was found to be destabilizing in the 
range of typical lead-lag and torsion frequencies. The main destabilizing 
factors were shown to be the bending-torsion structural coupling and, for 
configurations of low torsional stiffness, torsional dynamics. 
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Additional numerical results were presented by Friedmann in refer- 
ences 21-23. In references 21 and 22, flap-lag structural coupling was 
omitted and an error in the equations led to incorrect conclusions regarding 
the importance of blade droop (as mentioned in the discussion of ref. 22). In 
reference 23, flap-lag structural coupling was included and thus the problem 
treated was similar to that of reference 20, except for the representation of 
torsion motion. 

With the possible exception of reference 20, the work cited above does 
not furnish extensive numerical results for the stability characteristics of 
elastic cantilever rotor blades. This report is intended to provide results 
for a broad range of system parameters and to illustrate the general stability 
characteristics of hingeless rotors. It is an extension of reference 20 and 
incorporates the following changes: (1) the assumption of a blade model with 

a solid cross section is discarded, thereby permitting the stiffness and 
inertial parameters of the blade to be specified independently; (2) a suffi- 
ciently large number of assumed mode shapes is used to ensure convergence; 

(3) along with some notational changes, the blade configurations are defined 
in terms of rotating natural frequencies rather than nonrotating frequencies; 
and (4) blade configurations with various degrees of flap-lag structural 
coupling are represented by introducing a structural coupling parameter (R. 


Procedure 

The derivation of the general nonlinear equations of motion for an elastic 
rotating cantilever beam (ref. 1) is summarized in this report and specialized 
to a basic cantilever blade configuration without twist, nonuniformities in 
mass and stiffness, or chordwise offsets. A detailed derivation of the aero- 
dynamic forces is given based on strip theory and on a quasi-steady approxima- 
tion of unsteady aerodynamic theory for a two-dimensional airfoil. The non- 
linear partial differential equations are transformed into a system of ordinary 
differential equations by Galerkin’s method. These nonlinear equations of the 
generalized coordinates are linearized for small perturbation motions about 
the equilibrium operating condition. An iterative solution of the nonlinear 
algebraic equilibrium equations enables the coefficients of the perturbation 
differential equations to be determined. Finally, the eigenvalues of the per- 
turbation equations that determine the dynamic stability of blade motions 
about the equilibrium are obtained. Coupled mode shapes, determined from free 
vibrations {in vacuo) of the blade about the equilibrium operating condition, 
are used. 

The equations are solved for a variety of practical hingeless rotor con- 
figurations with emphasis on the variation of torsion frequency, lead-lag fre- 
quency, precone, the structural coupling parameter, and the collective pitch. 
Results are presented in the form of plots of equilibrium tip deflection, 
locus of roots, lead-lag damping, and stability boundaries. 
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GENERAL HINGELESS ROTOR CONFIGURATION PARAMETERS 


A comprehensive study of hingeless rotor stability is a formidable task 
because there is a large number of important configuration parameters. The 
purpose of the present report is to investigate the primary parameters of a 
somewhat simplified configuration in a systematic manner; the remaining param- 
eters are not considered. In order to place the present simplified configura- 
tion in proper perspective it is appropriate to describe the physical charac- 
teristics and define the parameters of a typical hingeless rotor blade. 

The general cantilever blade 
structure, shown in figure 1, is 
composed of two flexible beam seg- 
ments joined by the pitch change 
bearing. The inboard hub segment is 
fixed to the hub at the root end of 
the blade, while the outboard seg- 
ment can be rotated about the pitch 
change bearing by vertical movement 
of the pitch link from the swash- 
plate controls. Pitch link flexi- 
bility, represented by a spring ele- 
ment, will permit rigid body pitching 
motion of the outboard blade segment 
(i.e., root torsion). In general, 
both the inboard and outboard seg- 
ments of the rotor blade will have 
large nonuniformities of mass, bend- 
ing stiffness, and torsional rigidity, and both segments will bend and twist. 
These nonuniformities contribute to the important and complex structural 
coupling between the flap bending, lead-lag bending, and torsion of hingeless 
rotor blades. Some additional configuration parameters of importance are pre- 
cone, droop, torque offset, and sweep, all of which are illustrated in 
figure 2. These geometric parameters are usually tailored to minimize the 
steady state rotor blade stresses or to improve rotorcraft flying qualities, 
but they can also have an important influence on rotor blade stability. 

The above features may differ for certain special hingeless rotor con- 
figurations. For example, it is not uncommon for the inboard beam segment to 
be eliminated entirely by placing the pitch change bearing in the hub itself, 
thus leaving only a single outboard blade segment. Or the pitch bearing may 
be eliminated by making the inboard segment very flexible in torsion so that 
pitch change of the outboard segment is accommodated by twisting the inboard 
segment . 

The blade structure Itself may also influence the stability characteris- 
tics. A rotor blade may consist of several spar, skin, stiffener, and balance 
weight components, or it may be molded from many layers of composite materials 
bonded together. These materials must all be tailored to provide acceptable 
aerodynamic, strength, and stiffness properties. Additional parameters 


Rotor 

hub 



Figure 1.- Primary elements of a typical 
hingeless rotor blade (not all these 
elements are considered in the 
present analysis) . 
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influenced by these construction 
details are the positions of the mass 
center, tension center, elastic axis 
(shear center) , and the aerodynamic 
center. If these are not coincident, 
additional structural, inertial, and 
aerodynamic coupling between bending 
and torsion deflections will be intro- 
duced and may strongly influence rotor 
blade stability, as in the case of an 
articulated blade. The rotor blade 
may also be specially designed to 
produce desirable structural coupling 
effects; these effects may include 
equivalent pitch- flap and pitch-lag 
coupling, both of which can improve 
both rotor blade stability and rotor- 
craft flying qualities. Similar 
coupling effects may result from the 
kinematics of the pitch-link connec- 
tion between the outboard blade seg- 
ment and the swashplate control system. 
If the inboard blade segment is espe- 
cially long or flexible, this kine- 
matic coupling may be particularly 
significant because the motion of 
the pitch bearing due to blade bend- 
ing is increased. 



Droop 



Figure 2.- Configuration parameters 
of a hingeless rotor blade (of 
these parameters only precone is 
considered in the present analysis). 


As noted above, not all the rotor blade details described here are con- 
sidered in this report. The present analysis treats the elastic flap, lead- 
lag, and torsion deflections of a simple untwisted uniform blade with precone, 
but excludes pitch link flexibility, droop, sweep, or torque offset. Further- 
more, mass, tension axis, and aerodynamic center offsets from the elastic axis 
are all zero. The equations of motion are written for a single outboard blade 
segment without a flexible hub segment inboard of the pitch bearing. However, 
the effects of variable structural coupling due to the inboard and outboard 
segments are simulated by introducing the structural coupling parameter 
to arbitrarily vary the structural coupling in the uniform blade equations. 


EQUATIONS OF MOTION 


The equations of motion are taken from reference 1 and specialized for 
the present simplified blade configuration having uniform mass and stiffness, 
no twist, and no chordwise offsets of the elastic axis, tension axis, or 
center of mass. A brief outline of the derivation is given for completeness, 
and one modification to the equations is made to accommodate an approximate 
structural coupling parameter for expanding the scope of the numerical results. 
The expressions for aerodynamic forces are not included in reference 1; they 
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are derived in detail below on the basis of two-dimensional unsteady aerody- 
namic theory, 


Coordinate Systems 

Consider the undeformed beam in figure 3 with its elastic axis coincident 
with the X axis of the x, y, z coordinate system rotating with a constant 
angular velocity about a fixed point at the origin. The y axis lies in the 


Axis of 
rotation 



Figure 3.- Rotor blade coordinate sys- 
tems and deflections. 


z 



Figure 4,- Rotor blade cross section 
before and after deflections; note 
that y’ and z' are projections on 
the y 5 z plane. 


plane of rotation and the x axis 
is rotated through a small angle 
3pc from the plane of rotation. 
Bending deflections of the beam are 
defined by the displacements u, v, 
and w of the elastic axis parallel 
to the X, y, z coordinates, 
respectively, A second coordinate 
system, x’ , y’ , and z’, fixed to the 
blade with y', z’ axes parallel to 
the beam cross section principal 
axes, moves with the blade as it 
undergoes bending displacements , 
torsional displacements, and pitch 
angle rotation, 0, The projection 
of the beam cross section in the y, 
z plane is shown before and after 
deformation in figure 4. Before 
deformation, the blade principal axes 
are rotated with respect to the 
undeformed coordinates by the pitch 
angle. After deformation, the elas- 
tic axis is displaced by u, v, w 
and the blade twisted through the 
angle cf). Although the deformed 
y’ , z' axes do not lie exactly in 
the y, z plane, their projection 
in that plane is shown in figure 4. 

A small apparent twist, discussed in 
reference 24, is neglected by ignor- 
ing the difference between the actual 
and the projected twist angle (f) , 
but this must be accounted for in 
deriving the blade aerodynamic 
forces. 
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Structural and Inertial Terms 

The equations of motion are derived using Hamilton’s principle 


/.*^2 

j isu - 6T - 6W)dt = 0 
tl 


( 1 ) 


where the strain energy U, the kinetic energy Ty and the virtual work of 
external forces 6W are given by 


R 


^ i / // <Wxx + ^xn^xn + 

0 A 


° 1 / // pv • V dn d? dx 

0 A 


\ 


( 2 ) 


6W 


J* ^v^^ ^w^^ ^ M.6(j))dx 

0 ^ 


/ 


For the strain energy, the stresses are proportional to the strains; 

^xx ~ ^^xx> ^xri ” ^^xri’ ^xC “ ^^x^* kinetic energy, V is the 

velocity vector of an arbitrary point in the blade. The loads L^, L^, L^, 
and are of aerodynamic origin and are treated below. The strain displace-' 

ment relations as derived in reference 1 are: 


2 

Exx = + — + — + (e*4>’ + 2 




(n^ + - Xcf)” 


- v" [n cos(0 + (}))- ^ sin(0 + (j))] 

- w" [n sin(0 + (j)) + t cos(0 + (j))] 


^XTl 






(3) 


where X is the warp function for the cross section. The function A(n, C) 
may be determined by solving Laplace’s equation for the specific cross section. 


The remaining details of the derivation of the equations of motion are 
given in reference 1. An essential feature of the derivation is the introduc- 
tion of an ordering scheme in which e, a small parameter of the order of mag- 
nitude of the bending slopes, is introduced. The assumed order of magnitude 
of all quantities used in the equations appears in table 1. Within any equa- 
tion, terms of O(e^) are ordinarily dropped unless, for example, they 
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contribute to damping. The resulting dimensional equations are as follows 
where the radial aerodynamic force component is neglected as a higher 
order term. 


6u equation: -T’ - m^^^x - 2vaQ^7 = 0 

6v equation: 

-(Tv')' + [EI^. - (EI^. - ElyOsin^ 0]v"" + (EI^* - EIyt)[ ^^^” w"" 

- sin 26((j)v") ** + cos 26((|)w”) "J + 2mQ(u - 3pcw) + m(v - 
6w equation: 

-(Tw')' + [Ely. + (EI^. - EIy.)sin2 6]w"" + (EI^* - Ely . ) v"" 

+ cos 29((j)v")" + sin 29 (<j>w") " j + 2mfi6p(,v + mw = 

6(j) equation: 

(EI^. - EIy,)[^(w"2 - v"2) + v"w" cos 29 ]- GJ-ji" - kJCT^)')' 

+ ml4^ + mf22(k2^ - k2^)(|, cos 29 = “ *4^) ^^’^2 


TABLE 1.- ASSUMED ORDERS OF MAGNITUDE FOR PHYSICAL QUANTITIES 


II 



0(e) 

V w _ 

R ’ R ” 

0(e) 

a 



0(e2) 

II 

> |g 

0(e) 

d R 

J2dt ’ 

II 


0(1) 

^pc 

0(e) 

EA 



0(e-2) 

<t> = 

0(e) 

Ely, 

EIz' 

GJ 

= 0(1) 

0 = 

0(1) 


ii^Fr^’ 


*"a S 

R ’ R 

L 

’ R 

11 

0(e) 

0 (aerodynamic 
forces only) 

= 0(e) 




iiiil 1 II 

1 nil 1 1 «■■■ 

■■■III II 


(4a) 


(4b; 


(4c) 


(4d) 


mil 



where the tension is given by 


T = 


EA 


f' 



and 


(5) 


A = dn d? 
ly. E JJ ?2 dn de 

S 

= JJ* dn d? 


^v’ + ^z’ 

kA^ ^ 


A 


m E JJ p dn d? 
m 

^ m 

^2 = m JJ 


dn d? 




( 6 ) 


These cross-sectional integrals in equation (6) define the structural and mass 
properties of the rotor blade cross section. The subscripts s or m refer 
to integration over the portion of the cross section which contributes struc- 
tural or mass characteristics, respectively. These properties are the radial 
stiffness EA, the flapwise and chordwise bending stiffness Ely* and EI^’, 
the polar radius of gyration k^, the torsional rigidity GJ, the mass per 
unit length m, the flapwise and chordwise principal mass moments of inertia 
mk^^ and mk^^, and the torsional mass moment of inertia nik^. 

The present nonlinear equations may be compared directly with the well- 
known linear equations of Houbolt and Brooks (ref. 2). When specialized to 
untwisted beams with uniform mass and stiffness and with no chordwise offsets, 
the equations of Houbolt and Brooks are identical to the present equations 
with the exception that the underlined terms in equations (4) are not present. 
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The underlined nonlinear terms arise from both inertial and structural effects 
and they are important for dynamic analyses of cantilever rotor blades. (The 
underlined terms due to precone 3^^ are linear, but precone was not consid- 
ered by Houbolt and Brooks.) 

The various terms in the present equations will now be briefly described. 
We first discuss the inertial terms, and in particular the nonlinear flap-lag 
inertial coupling terms. In the 6u equation, the term -2mS7v represents 
the centrifugal force of lead-lag deflection velocity. This dependence of 
tension on lead-lag deflection produces nonlinear products when T is substi- 
tuted into the (Tv’)’ and (Tw’)’ terms in the lead-lag and flap bending equa- 
tions. A similar nonlinear expression for the 2m^^u Coriolis term in the 
lead-lag equation may be obtained by solving for u in terms of v, w, and T 
in the nonlinear strain relation, equation (5), As discussed in reference 6, 
these terms must be included to properly represent the Coriolis and centrifugal 
coupling between the flap and lead-lag equations. In the torsion equation, the 
tension-torsion coupling term (Tcj) ’ ) ’ , has a small coefficient and only the 

linear contribution -k^^ (m^^^/2) [ (R^ - x^)(j)’]’ is retained. Inertia terms in 
the torsion equation include the torsional inertia mk^^$ and the tennis 
racket effect; the tennis racket effect is proportional to the difference 
between the principal mass moments of inertia. We next discuss the structural 
terms. 

In addition to the familiar torsional rigidity term -GJcf)", there are the 
nonlinear bending-torsion structural coupling terms underlined in the torsion 
equation, (eq. (4d)). These terms are probably the most important ones for 
hingeless rotor stability because they are large and sensitive to many config- 
uration parameters. They are actually moments proportional to the difference 
between lead-lag and flap bending stiffness EI 2 ' " Elyf, and products of the 
flap and lead-lag curvatures. These terms are discussed by Mil’ (ref. 9) as 
noted in the introduction; the present equations also include the correspond- 
ing nonlinear bending-torsion terms in the flap and lead-lag bending equations 
(eqs. (4b), (4c)). The flap-lag structural coupling terms in equations (4b) 
and (4c), with coefficients proportional to EI^’ - EIy» are linear and are 
included in the equations of Houbolt and Brooks. These terms couple flap bend- 
ing and lead-lag bending elastically when 0 0 and, like the bending-torsion 

coupling terms, they are proportional to the difference between the lead-lag 
and flap bending stiffnesses. In some analyses, such as those of references 7, 
18, 21, and 22, these terms have been discarded, although they strongly influ- 
ence the stability characteristics of some hingeless rotor blade configurations. 

In general, the flap-lag structural coupling terms are also dependent on 
blade twist and on bending stiffness distributions in the case of nonuniform 
blade configurations. In particular, the flap-lag structural coupling is 
dependent on the relative stiffness of the blade segments inboard and outboard 
of the pitch bearing. This is because the principal elastic axes of the out- 
board blade segment rotate through an angle 0 as the blade pitch varies 
while the inboard segment principal axes do not- The resultant effective 
orientation of principal axes depends on the blade geometry and distribution 
of bending stiffness inboard and outboard of the pitch bearing. Although the 
variations in these coupling effects significantly influence stability, they 
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are not present in the case of a simple single segment uniform beam; they are 
difficult to include exactly without resorting to a more general blade con- 
figuration and a more sophisticated analysis- However, an approximate repre- 
sentation of these effects may be introduced in the present equations with no 
increase in complexity. This is accomplished by arbitrarily assuming that the 
average inclinatipn of the principal elastic axes of a nonuniform blade is 
equal to some fraction of the inclination of the principal axes of a uniform 
single segment blade. This entails replacing 0 by (R0 in the structural 
terms in the equations while the mass and inertial terms are unchanged- The 
factor (R is called the structural coupling parameter. When (R - 1 , the 
original equations are retained, but as (R is reduced to zero, the flap-lag 
structural coupling terms diminish and eventually vanish. The structural cou- 
pling parameter (R also influences the bending-torsion structural coupling 
terms , but to a lesser degree . 

As defined here, the structural coupling parameter (R is nearly equiva- 
lent (for small pitch angles) to an analogous parameter R used in refer- 
ence 5 to characterize the variable flap-lag structural coupling of the 
approximate rigid, hinged blade representation of an elastic blade. Additional 
discussion of flap-lag structural coupling is given in reference 5. 

In their present form, including the nonlinear strain relation, there are 
five nonlinear equations in u, v, w, c|), and T. These equations can be sim- 
plified by substitution to eliminate u and T as dependent variables, follow- 
ing the procedure given in references 1, 6, 19, and 20. First, the 6u equa- 
tion is integrated to obtain the tension: 

f)2 ^ 

T = (R^ - x^) + 2mQ J V dx (7) 

X 

This equation is substituted into the 6v, 6w, and dcj) equations to eliminate 

R 

T. When equation (7) is substituted into the dcf) equation, the 2m^^ J* v dx 

X 

term is discarded because it is O(e^). The 2mQu term in the 6v equation 
is replaced as follows, solving first for u’ from equation (5) and substitut- 
ing equation (7) for T. 

, m^2^R^ fl - x^/R^\ , 2m^^^R^ v dx v'^ w*^ 

“ ■ V — 1 J 155 — 2 r 

This is integrated with respect to x and differentiated with respect to t 
yielding 


EA 






L 

1 dx dx - 

__ n o 

J (v’v* + w’w’)dx 

(9) 


0 ^ 0 


While not essential, it is convenient to discard the first term of equation (9) 
since mJ2^R^/EA = O(e^) where e is of the order of magnitude of the bending 
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slopes. This is equivalent to postulating that the rotor blade is inexten- 
sional for perturbation bending deflections; that is, radial shortening of the 
beam is a purely geometric consequence of the transverse bending deflections 
of the blade. This assumption is sometimes made on the basis of intuitive 
reasoning; however, it is strictly valid only for sufficiently large values of 
the dimensionless radial stiffness, EA/mf2^R^. 

After substituting equation (9) into the 6v equation (4b) , three equa- 
tions for V, w, and (j) remain. In the structural terms, 0 has been replaced 
by ^0 according to the discussion above. 


6v equation: 


+ 



sin(2(R0) 

(EI^* ” ^^y ^ j ^ 


- 2m^^3p^w - 2 m^^ J (v 


+ [EI^* - (EI^' - EIy.)sin2(«e)]v"" 
y,)[-s±n(2iR0)(((>v”)” + cos(2(R0)(4>w”)”] 
'v' + w'w')dx + m(v - = L-^ (10a) 


6w equation: 


[w’ (R^ - x^)]’ - 2m^2 


J V dx^ 


+ [Ely. + (EI^.' - EIy.)sin2((R0)]w"" 


+ (EI^, - Ely.) v"" + (EI^. - EIy.)[cos(2(R0)(<j>v")" 


+ sin(2(R0) (^w”)”] + 2m^^3p^v + mw = (10b) 


pc^ 


6(j) equation: 

- ^ k/[.|>’(R2-x2)]’ - GJ(|)" + (El^. - Ely.) j^(w"2 

+ v"w" cos(2(R0)l + ~ cos 20 


„2x sin(2(R0) 

V ) ^ 


= - 14 ^) 


s in ( 2 0 ) 


(10c) 


These equations are nonlinear, integro-partial differential equations with 
variable coefficients in x. In the next section, expressions for L^, L^, 
and M(|j, the aerodynamic loading terms, will be derived. 
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Aerodynamic Loading 


The aerodynamic lift and pitching moment acting on the blade in hover are 
based on Greenberg’s extension of Theodorsen’s theory (ref. 25) for a two- 
dimensional airfoil undergoing sinusoidal motion in pulsating incompressible 
flow. The rotor blade aerodynamic forces are formulated from strip theory in 
which only the velocity component perpendicular to the blade spanwise axis 
(the x’-axis in the deformed blade coordinate system x’ , y’ , z’ in fig. 5) 
influences the aerodynamic forces. A quasi-steady approximation of the 
unsteady theory for low reduced frequency k is employed in which the 
Theodorsen function C (k) is taken to be unity. The steady induced inflow for 
the rotor is calculated from classical blade element-momentum theory. These 
simplifying assumptions are judged to be adequate for low frequency (mainly 
determined by the blade bending frequencies) stability analyses of a hovering 
rotor. 


In Theodorsen’s theory (ref. 26), 
a two-dimensional airfoil is assumed 
to be pivoted about an axis which may 
be distinct, in general, from the 
aerodynamic center axis. The airfoil 
is pitched at an angle s(t) to the 
free stream flowing at constant 
velocity V. The airfoil is verti- 
cally displaced with velocity h(t) 
positive downward as shown in fig- 
ure 5. Greenberg has extended 
Theodorsen’s theory for pulsating 
free- stream velocity V(t). The 
relations for lift and pitching 
moment per unit length may be 
expressed in terms of the circulatory 
and noncirculatory components 



Figure 5.- Rotor blade airfoil section 
in general unsteady motion. 
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Sc 


= S + Sc 


( 11 ) 


With the airfoil pivot axis (analogous to the rotor blade elastic axis) at the 
airfoil quarter chord (the airfoil aerodynamic center) these components are 
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NC 
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p ac 
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p ac 
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+ Ve + -| 
Poo'^^/c\3 e 

2 Uj 2 


p ac 

~r~ 



2 

Ve 


( 12 ) 
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The Theodorsen function C(k) has been set equal to unity in the circulatory 
lift. It should be noted that e is the angular position of the airfoil with 
respect to space; e and e are the angular velocity and angular acceleration 
of the airfoil. The instantaneous angle of attack of the airfoil 
a = tan“^ (Up/U-p) is the angle between the airfoil chord line and the resultant 
fluid velocity U of the airfoil. The airfoil velocity components in the 
principal axis system y’ , z’ are U-p and Up shoxm in figure 5. It is 
desirable to express the aerodynamic forces and moments in terms of Up and Up. 
Assuming that the angles e and a are small yields 


Up = - h - Ve 


u = /u-r^ + Up2 s V 


I (13) 


Substitution of equations (13) into equations (12) yields 


p ac 


NC 


-f (- 


Uo + 


p ac 




(- 


U -Ut3 + 


f^) 


(14) 


Next we consider the total aerodynamic forces in directions parallel and 
perpendicular to the airfoil chord line. The noncirculatory lift is taken to 
act normal to the chordline, and the circulatory lift is taken to act normal 
to the resultant blade velocity U. An aerodynamic profile drag force per 
unit length, acting parallel to the resultant blade velocity, is included 
based on a constant profile drag coefficient c^^ . 


p ac 


D = 


'0 


(U^^ + Up^) 


(15) 


The force components and directions are shown in figure 6. The force compo- 
nents T, normal to the airfoil chord line, and S, parallel to the airfoil 
chord line , are therefore 


From figure 5, 


T = cos a + + D sin a 


S = -Lq sin a - D cos a 


cos a = 


Ut 


u 


+ Up2 


Up Up 


sin a = 


U 


+ Up2 


(16) 
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Figure 6.- Orientation of components of aerod 3 mamic loading. 

Substitution of equations (14), (15), and (17) into equations (16), with 

CA /a neglected with respect to unity, yields 
0 



p ac 
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The expressions for aerodynamic pitching moment components may be written from 
equations (12) and (13) as 


p ac 

GO 


M 


c 


p ac 

oo 

2 ~ 



+ 


3c 




(19) 


/ 


where U has been approximated by U-j- in Mq. The total pitching moment is 
then given by 


M 


<!> 


M = 


p ac 

cx> 

2 



( 20 ) 


In appendix A, Up, Up, and t are expressed in terms of 
the appropriate order of magnitude. From equations (A3) and 


V, w, and (j) 
(A9) we have 


to 
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Up = -^^xle + (f) + 


J" v’w" dxj - (0 + (|))v + + w + ^v(3p^ + w’) \ 


Up = f^x 4- V 
e = 3> + S^(3p(, + w') 


> ( 21 ) 


where is the induced inflow velocity. The underlined terms in Up pro- 

duce second order terms in the angle of attack; they are discussed in detail in 
reference 24. These terms have been shown to influence significantly the flap- 
lag-torsion stability of cantilever rotor blades (ref. 20). The transformation 
[T] may also be used to resolve the blade forces S and T into and 

parallel to the y and z axes of the undeformed blade coordinate system. 
Retaining terms to the appropriate order of magnitude yields 


= S - T(0 + cf)) 
= T + S(0 + (|)) 


(22) 


Substitution of equations (21) into equation (20) and equations (18) and (21) 
into equations (22) yields the following expressions for L^, L^, and M: 


p ac 
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^ - Q^x^ ^ - $7xvp(0 + <j)) - 


^0 
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+ [2vp - Qx(0 + (j))]w 


2f2x + (0 
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+ <j))v^ V I 


oo 
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^0 + ({) + J* v'w"^ - fi^xv(Bp^ + w') 


+ (®pc + [2nx(0 + (f>) - v^]v - Qxw 


3c • c 
+ ^ nx<|) - w 


"I 


p ac 




(x. 


(23) 


Nonlinear rate product terms w, v^, and w^ are neglected since they do not 
contribute in a linearized stability analysis. Also, all O(e^) terms, except 
those that contribute to lead-lag or torsion damping (double underlined) , are 
neglected (including e terms). The single underlined terms are the second 
order angle of attack terms noted above. The induced inflow v^ is taken to 
be steady and uniform along the blade radius equal to the value of nonuniform 
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inflow given by blade element momentum theory at the radial station x = 0.75R 
(ref. 27). The blade angle at x = 0.75R is set equal to the blade collec- 
tive pitch plus the equilibrium elastic twist <()q at x = 0.75R. Thus 


Vj_ = sgn[0 + (|)q(0.75R)]^^R 


8 y * iro 


e + <|)„(0.75R) 




(24) 


where a is the blade solidity bc/irR. 


Solution of the Equations of Motion 


In this section we describe the solution of the final flap-lag-torsion 
equations of motion obtained by combining equations (10) and (23). These non- 
linear, variable coefficient, integro-partial differential equations are 
solved by Galerkin’s method using six coupled, rotating mode shapes to evalu- 
ate the stability of small perturbation motions about the equilibrium operating 
condition. These coupled, rotating mode shapes are determined by first calcu- 
lating the equilibrium position and formulating the stability analysis in 
terms of standard, uncoupled, nonrotating cantilever beam mode shapes. The 
coupled mode shapes are then determined from an approximate free vibration 
analysis (in vacuo) of the rotating beam about the deformed shape of the equi- 
librium operating condition. 


We begin by first reducing the equations of motion to ordinary differen- 
tial equations by Galerkin’s method. The dimensionless bending (v = v/R, 
w = w/R) and torsion deflections are first expressed in terms of a series of 
generalized coordinates and mode shape functions: 


N 

V = Yi V. 

j=l 

N 

w = 2 (<1')'*'^ (x) 

j = l 


\ 


> (25) 


N 

<j> = 23 (^) 

j = l 

where ij; = ^^t and x = x/R. This operation yields 3N modal equations in 
terms of modal generalized coordinates Vj , Wj , and Oj (appendix B contains 
the results of this operation) . The assumed mode shapes for the bending and 
torsion deflections are the standard nonrotating, uncoupled mode shapes for a 
uniform cantilever beam: 


(x) = cosh(3jx) - cos(3jx) - aj[sinh(3jx) - sin(3jx)] 

04 (x) == /l sin(y.x) 

J 


i (26) 
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The constants aj and 3j are tabulated in reference 28 and Yj = *^(3 ■“ 1/2). 
Nonrotating mode shapes are used because of computational ease. Since they 
depend only on the fixed constants aj and 3j > and not on EIy» and the 

modal integrals that result when Galerkin’s method is applied need to be cal- 
culated only once (see refs. 6 and 19, for example). Furthermore, many of the 
modal integrals have closed form solutions in terms of aj and 3j (ref. 29). 


The 3N nonlinear, nonhomogeneous , constant-coefficient ordinary dif- 
ferential equations are then linearized for small perturbation motions about 
the equilibrium operating condition by expressing the time-dependent general- 
ized coordinates in terms of steady equilibrium quantities and small unsteady 
perturbation quantities. 

Vj w = Vqj + AVj W ) 


Wj W = + AWj W 


(27) 


j 


Two sets of equations are obtained from this operation. First, substituting 
only the steady equilibrium quantities into the nonlinear modal equations 
yields 3N nonlinear algebraic equations in Vqj , ^0 j > (given in 

appendix B) which define the equilibrium deflections. These equilibrium equa- 
tions are solved by the Newton-Raphson method. The second set of 3N 
equations is obtained by substituting equations (27) into the modal equations, 
subtracting the equilibrium equations, and discarding all nonlinear products 
of perturbation quantities. The coefficients of these linearized perturbation 
equations are functions of the equilibrium solution. The perturbation equa- 
tions define the unsteady blade motion near the equilibrium operating condition 
and the stability of this motion is determined using standard techniques. The 
perturbation equations are linear, homogeneous, constant-coefficient ordinary 
differential equations of the form 


[M]{X} + [C]{X} + [K]{X} = 0 


(28) 


where [M] is symmetric, and where [C] and [K] are asymmetric and depend on 
Voj ,W q j , and ^oj • Each matrix is 3N x 3N (these matrices are given in appen- 
dix B)'^and the vector {X} is of length 3N and is given by 


I AVj ’ 
AWj 
A$. 


(29) 


The stability of motion about the equilibrium operating condition is deter- 
mined by the eigenvalues of the 6N x 6N matrix [P] where 
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X 




-M~lc 



= [P] 


X 

1 

-M"^K 


X 

1 1 

X 


Since we are primarily concerned with lower frequency instabilities (first 
lead- lag, first flap, and first torsion frequencies), there is a value of N 
for which any increase in N will not appreciably change the eigenvalues 
associated with these lower frequencies. It is at this value of N that the 
eigenvalues are considered to be converged. For practical hingeless rotor 
configurations, N = 5 gives suitably converged eigenvalues; the matrix [P] 
is thus 30 X 30. By a change of modal coordinates, the size of the matrix [P] 
may be greatly reduced without significantly changing the eigenvalues of 
interest. Such a transformation may be found by first considering free vibra- 
tions (in vaouo) of the blade about the equilibrium deflected state. The 
equation of motion, analogous to equation (28), is 

[M]{X} + [G]{X} + [K^]{X} - 0 (31) 

where the subscript v implies the vacuum case. Both [M] and 
symmetric; [K^] is equal to [K] with all aerodynamic terms set equal to 
zero, and [G] is antisymmetric and equal to [C] with all aerodynamic terms 
set equal to zero. The presence of the matrix [G] causes the eigenvectors 
of the free vibration to be complex. This may be avoided for computational 
efficiency by approximating equation (31) as 

[M]{X} + [K^]{X} = 0 (32) 


The matrix of eigenvectors [U] for equation (32) is orthogonal with respect 
to [M]. Thus, 

[U]^[M][U] = [I] (33) 


According to Meirovitch (ref. 30) a so-called principal coordinate transforma- 
tion for equation (28) may be determined by replacing {X} by [U]{X}. We 
may then pre-multiply equation (28) by [U]"^ to take advantage of the form of 
equation (33) yielding 

[I]{X> + [u'^CUHX} + [u\u]{X} = 0 (34) 


Hence , 
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0 

I 

( X 

1 

X 

X ' 


_-u%u 

-u'^cu. 

(X 

= [p*] 

X 


(35) 


The matrices [P] and [P*] have the same eigenvalues. However, because of 
the nature of this modal coordinate transformation from [P] to [P*] , the rows 
and columns corresponding to high frequency modes of both [u'^KU] and [U*^CU] 
may be removed without affecting the eigenvalues of the low frequency modes of 
interest. These 3N x 3N matrices are thus reduced to M x M matrices whose 
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rows and columns correspond to the M low frequency modes that are retained. 
The rows and columns that are retained in [u'^KU] and [u'^CU] may be chosen in 
two ways: (1) the M rows and columns that correspond to the M lowest fre- 

quency modes of the blade may be retained, or (2) the M rows and columns that 
correspond to M modes selected arbitrarily from the lowest lead-lag, the 
lowest flap, and the lowest torsion frequency modes are retained. For the 
second case, under certain conditions, M = 3, 4, or 5 will result in converged 
eigenvalues. In either case, suitably converged results do not require M > 6 

The reduced matrices are analogous to stiffness and damping matrices 
generated from M coupled, rotating modes. Since the analysis is formulated 
in terms of standard cantilever mode shapes, however, repeated numerical 
integration of modal integrals is not necessary for different values of beam 
stiffnesses. Instead, the matrix operations described above lead to a 
net savings in CPU time. For all results given in this report, M = 6 unless 
otherwise noted. Thus [P*] is reduced from 6N x 6N (30 x 30) to 2M x 2M 
(12 X 12). Eigenvalues may be calculated from a 12 x 12 matrix in less than 
10 percent of the time required for a 30 x 30. The total savings in CPU time 
is about 60 percent for each case including the equilibrium solution, the 
above modal analysis, and the eigenvalues for the perturbation equations. 


Approximate Equations of Motion 

The flap-lag-torsion equations may be considered as an extension of the 
earlier flap-lag equations of reference 6. They permit the analysis of tor- 
sionally flexible rotor blades, but also introduce additional equations and 
degrees of freedom for the torsion generalized coordinates. It is possible to 
reduce the flap-lag-torsion equations to an approximate set of modified 
flap- lag equations and still retain the most important effects of torsional 
flexibility. This is possible because the structural terms in the torsion 
equation are usually large compared to the inertial and aerodynamic damping 
terms. Based on order of magnitude considerations, the and Oj terms are 

small in comparison to the Oj terms when is sufficiently large. The 

approximation simply consists of discarding the $j and terms in the 

torsion modal equations leaving (in addition to the tension-torsion and tennis 
racket terms) the torsion rigidity term proportional to and the bending- 

torsion structural terms proportional to products of the bending generalized 
coordinates Vj and Wj . This permits to be expressed in terms of Vj and 

Wj and permits Oj to be substituted into the flap and lead-lag modal equa- 
tions. This reduces the number of equations and eliminates the torsion gener- 
alized coordinates Oj as independent degrees of freedom. The accuracy and 
limits of validity of the resulting modified flap-lag equations will be 
examined below. 


SELECTION OF CONFIGURATION PARAMETERS 


A few brief comments are in order to explain clearly the choice of the 
blade configuration parameters. These consist of the structural and inertia 
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parameters Ely t , EI^ » , GJ/m^2^R^, Ra/^h* ^m;^/^m 2 > and <R, and 

the geometric and aerodynamic parameters c/R, 3pc» and a, and the Lock number 
y = 3p^acR/m. The pitch angle 6 is considered to be the control parameter 
defining the rotor operating condition; that is, the rotor thrust. 

Beginning with the geometric and aerodynamic parameters, typical values 
are chosen as follows: the Lock number y = 5.0 and the solidity a = 0.10. 

Since a - Lc/ttR, the choice of a four blade rotor (b = 4) yields the dimen- 
sionless chord ratio c/R = tt/ 40. The geometric blade precone angle gpc has 
an important influence on stability, and thus a range of values will be 
investigated. 


The structural and inertial parameters mainly determine the natural fre- 
quencies of the blade which in turn serve as rotor blade configuration param- 
eters. In general, explicit expressions do not exist for the structural and 
inertial parameters in terms of the blade natural frequencies, and it is 
necessary to solve for these parameters iteratively if the natural frequencies 
are specified. In this report, the following procedure is used. For 
an untwisted cantilever beam with uniform mass and stiffness, the structural 
and inertial parameters are first expressed explicitly in terms of the uncou- 
pled natural frequencies for the nonrotating condition. These relationships 
are 


EIz' 


V. 


NR 


3^4 ’ 


Ely. 

mfi2R4 




W, 


NR 


GJ 

mf^^k^R^ 


^NR 


(36) 


where the subscript NR refers to 
frequencies of the nonrotating beam. 
The nonrotating frequencies are then 
used as an intermediate step in deter- 
mining the uncoupled (at 0=0), 
rotating natural frequencies. For an 
untwisted uniform cantilever beam, the 
rotating bending frequencies may be 
accurately expressed in terms of the 
nonrotating frequencies using the 
approximate analytical relations 
developed in reference 31. These 
relationships are shown in figures 7 
and 8. An analogous relationship for 
the rotating and nonrotating torsion 
frequencies, also derived in refer- 
ence 31, depends on the intensity of 
tension-torsion coupling that is 
proportional to 

inertia ratio This result 

is shown in figure 9 based on a 
choice of 

^mi/^m 2 ™ 0.0. The remaining torsion 
inertia parameter, the mass radius of 



Figure 7.- Rotating lead-lag fre- 
quency versus nonrotating lead-lag 
frequency: 6 = 0.0. 
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Figure 8.- Rotating flap frequency 
versus nonrotating flap frequency: 

e = 0. 



Figure 9.- Rotating torsion frequency 
versus nonrotating torsion frequency; 
9 = O.Oj ~ 0.0, and 

~ 1.5. 


gyration k^/R, is taken to be 0.025. 
Therefore, any choice of the three 
fundamental rotating frequencies, 
together with the values for 

given 

here and the frequency variations of 
figures 7-9, will yield the appro- 
priate values of EIy»/m^^R^, 

EI2 ' , GJ/mQ^R^, kA^/km^, 
k^jj^/kjj^, and 1%/R for the equations 
of motion. 

Values for the blade natural 
frequencies will be chosen as fol- 
lows. The dimensionless flap natural 
frequency typically lies between 1.1 
and 1.2; for the present results 

is taken to be 1.15. The dimen- 
sionless lead-lag natural frequency 
^ typically falls into one of two 
ranges: 0.5 to 0.75 for soft inplane 

(o)^ < 1.0) configurations and 1.2 to 
1.7 for stiff inplane (o)^ > 1.0) con- 
figurations. Two typical values, 

03^ = 0.7 and 1.5, are chosen here. 
Since Elyi/EI^* = 0(1) was assumed 
from the outset in deriving the 
equations, the validity of the equa- 
tions is limited to < 3 (i.e., 

the practical range) . To analyze 
configurations with larger one 

should refer to reference 1 and 
rederive the equations based on a 
different ordering scheme, such as 
= 0(e“^). This will 
introduce additional terms and the 
equations will then be valid for 
larger values of oj^. The dimension- 
less torsion natural frequency 
usually ranges from 2.5 to 8.0. We 
will use 5.0 as a typical value and 
include results for other values as 
well. 


The structural coupling parameter (R will arbitrarily be given various 
values between 0.0 and 1.0. It is not possible to relate this parameter in a 
simple way to the stiffness distribution characteristics of actual rotor 
blades and the results can only indicate the range of behavior that may be 
encountered. 

In solving the equations using assumed mode shapes, the modal parameters 
are important. Unless otherwise specified, six coupled rotating mode shapes 
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(M = 6) will be used, thus producing six degrees of freedom. These are based 
on the free vibration analysis using 15 (3N) uncoupled nonrotating mode shapes. 
The six coupled rotating mode shapes retained are the six lowest frequency 
modes. The various parameter values discussed above are summarized in table 2. 

TABLE 2.- VALUES OF CONFIGURATION AND OPERATING CONDITION PARAMETERS 

USED FOR NUMERICAL RESULTS 



Baseline values 

Additional values 

0)v 

0.7, 1-5 

Variable 

% 

1. 15 

— 


2.5, 5.0, 8.0 

Variable 


1.5 

— 

0 , rad 

o 

o 

Variable 0.0 - 0.5 

3pc> rad 

0.0, 0.05, 0.1 

Variable 

(R 

0.0, 0.4, 1.0 

Variable 0.0 - 1.0 

vr 

0.025 

— 


0.0 

— 

Y = (3p^acR/m) ! 

5 

— 

o 

0. 10 

— 

c/R 

tt/ 40 

— 

b 

4 

— 


0.01/2TT 

— 

M 

6 

2-5 


RESULTS 


Steady State Deflections 

The steady state equilibrium deflections of the rotor blade must be deter- 
mined because they contribute to the coefficients of the linearized perturba- 
tion equations. The equilibrium deflections are required also to define pre- 
cisely the operating condition of the rotor. This is because the total blade 
angle, which determines the rotor thrust in hover, is equal to the blade pitch 
angle 0 plus the equilibrium elastic torsional deflection of the blade ())q. 

A few results are given to indicate the general behavior of the equilibrium 
deflections and the influence of several configuration parameters. The results 
are given in terms of the dimensionless tip deflections that are related to 
the equilibrium generalized coordinates by the following relations: 
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( 37 ) 


The blade pitch angle 0 is varied throughout the nominal range of thrust 
including very high thrust conditions (0 > 0.3 implies that C-p/a > 0.17). 

The results given in figures 10 and 11 show the equilibrium deflections 
of a soft and a stiff inplane rotor blade as a function of the pitch angle 9 
for several values of the structural coupling parameter (R. Increasing 6 
increases the aerodynamic forces acting on the blade and generally increases 
the magnitude of the equilibrium deflections. The structural coupling param- 
eter (R has little direct effect on the equilibrium flap deflection Wp 


tip’ 


although the variation of ^o|tip change the blade angle and thus 

indirectly influences ^o|tip* lead lag deflection is directly influenced 

by ^ as a result of the flap-lag structural coupling. When (R 0 and 

EI 2 ’ > ^^y’ > the blade tends to bend in a direction normal to the airfoil 

3 j_1 ^ j_ J — I j- ^ 3 3 T ; J — I 

I tip 


chord so that negative 

deflections when 0 > 0. 
ingly negative values of 


deflections are produced by positive w 


tip 


This is evident in figures 10 and 11 by the increas- 


tip 


as (R varies from 0 to 1.0. 


The equilibrium torsion deflection cj) 


tip 


depends on a combination of 


inertial and structural moments. First the inertial term, or tennis racket 
effect, that is independent of (R dominates the torsion moments at small pitch 
angles. This term produces a negative equilibrium torsion deflection propor- 
tional to 0 at small pitch angles that is clearly evident in figures 10 
and 11. Second, the structural torsion moments strongly influence the equi- 
librium torsion deflection, but in a complex manner. In the torsion equation, 
the bending-torsion structural coupling terms are nonlinear products of the 
bending deflections. Thus, parameters that influence the equilibrium values 


of 


'0 


tip 


and w. 


tip 


similarly influence the equilibrium torsion 
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(a) Flap bending. (a) Flap bending. 

(b) Lead-lag bending. (b) Lead-lag bending. 

(c) Torsion. (c) Torsion. 

Figure 10.- The effect of pitch angle Figure 11.- The effect of pitch angle 
and (R on the equilibrium deflec- and <R on the equilibrium deflec- 
tions of a soft inplane rotor blade: tlons of a stiff inplane rotor 

(ji = 0.7, 0 ), = 5.0, 3 = 0.0 rad. blade: to = 1.5, to = 5.0, 

^ 3 = O.O^rad. * 


deflection. To help explain the structural coupling terms that influence the 
equilibrium torsional deflection, the structural terms of equation (10) are 
written below. 


[EI^. 


(EI^, - Ely 

+ (EI^. 


)sin2((R0)]v"" + (EI^' - 
- Ely.) [-sin(2(Re)(<i)v")" 


■pi-r \ sln(2(R0) itii 

Ely.) 2 ^ 

+ cos (2<R0) ((fiw") "] + . . . 


(EIz- 


-GJcj)" 


- Ely.) v"" + [Ely. + (EIz> - EIy.)sln2((R0)]w"" 

+ (EI^' - Ely.) [cos(2(R0) (c()v")" + sin(2(R0)((J)w")"] + . . . 

+ (EI^. - Ely.) r(w"2 - v"2) + v"w" cos(2(R0) ~| + . . . 


\ 


> (38) 


Consider first the case where (R = 0 and the flap-lag structural coupling 
terms vanish. In the torsion equation, the bending-torsion structural moment 
(singly underlined in eq, (38)) is proportional to the product of the flap and 
lead-lag bending curvatures. This is a positive torsion moment that counter- 
acts the equilibrium torsion deflection due to the inertial tennis racket 
effect shown in figures 10 and 11. Furthermore, since the bending-torsion 
structural coupling is proportional to - EIy» this effect is larger for 

the stiff-inplane than for the soft-inplane configuration. In fact, for suf- 
ficiently large 0 or the structural torsion moments can produce a static 

divergence for low values of R. For the case when (R ^ 0, there is another 
torsional moment (doubly underlined in eq. (38)) that is proportional to 
squares of the bending curvatures. This moment generally acts to reduce the 
effect of its companion term; this is evident in figures 10 and 11 as R 
varies from 0 to 1. 


The two terms together are equivalent to a single product term 
(EI^* “ EIy»)Vp"wp'* where Vp and Wp are blade bending deflections in the 
direction of the blade principal axes. This axis system is obtained simply by 
rotating the y,z axis system through an angle R0. Alternatively, the non- 
linear structural torsion moment can also be shown to be proportional to the 
angle between the resultant direction of the blade bending deflection and the 
principal elastic axis angle R0. As (R increases, this angle, the positive 

structural torsion moment, and reduced. 

^ ' txp 


The effects of precone on the equilibrium deflections for a soft inplane 
and a stiff inplane rotor blade are shown in figures 12 and 13, respectively. 
Positive precone reduces the elastic flap deflection of the blade because of 
centrifugal force. For this example, R = 1.0 and flap-lag structural coupling 
with positive precone generates a positive lead-lag deflection increment due 
to the negative increment of elastic flap deflection. Similarly, the 
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0, rad 


(a) Flap-bending. (a) Flap-bending. 

(b) Lead-lag bending. (b) Lead-lag bending. 

(c) Torsion. (c) Torsion. 


Figure 12.- The effect of pitch angle 
and precone on the equilibrium 
deflections of a soft inplane rotor 


blade: 


= 0.7, 0). = 5.0, (R = 1.0. 


Figure 13.- The effect of pitch angle 
and precone on the equilibrium 
deflections of a stiff inplane 
rotor blade: o) = 1.5, = 5.0, 

«= 1 . 0 . ^ 


structural torsion moments, and hence the equilibrium torsion deflections, is 
influenced by precone. 


As noted above, the nonlinear equilibrium equations were solved by the 
Newton-Raphson method. In some cases, multiple equilibrium solutions were 
found, but in all cases only the physically relevant solution was statically 
stable. At high pitch angles, care was required to insure that the numerical 
scheme converged to the proper solution. This was done by using a previous 
solution at a lower pitch angle as the initial condition for a solution at a 
slightly higher pitch angle. 


Mechanism of Bending-Torsion Structural Coupling 
The equilibrium results showed that the torsion deflection 


tip 


xs 


dependent on ^^^^tip nonlinear structural bending- 

torsion coupling. This coupling also influences the linearized perturbation 
equations, producing perturbation torsion deflections proportional to pertur- 
bation flap and lead-lag bending deflections. This type of coupling is known 
(ref. 5) to have a strong influence on hingeless rotor blade stability and, 
therefore, a brief discussion of the equivalent kinematic pitch-flap and pitch- 
lag couplings generated by the structural coupling is included here. 


For simplicity, we consider the perturbation torsion equation for one 
flap, one lead- lag, and one torsion mode (N = 1) , neglecting torsional inertia 
and damping: 


- k^ 
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where 




m. 


CD.^ = 


V 




N 


GJtt" 




and the modal integrals defined in appendix B are 

1 


K 


111 


N 


11 


/ dx s 5.03911 

0 




2 dx = 1.07247 


30 


r nil 


I III IN I 


III I II null ■■ II I 


III I II I 


I nil I 



Taking the ratios of torsion deflection to flap and lead-lag bending deflec- 
tions yields, for =0, 


0O - 




^ 3A4>j 




3AW, 


-KiiiCEIz' - Ely.) 

(o)(j,^ - 1 + cos 20) 

^ ^ ^ -Kiii( EI^. - Ely.) 
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[Vqt cos(2(R0) + Wqi sln(2(R0)] 


[-Vqj sin(2(R0) + cos(2<R0)] 


(40) 


where 0g and 0^ are defined as the equivalent kinematic pitch-lag and pitch- 


flap couplings due to the structural properties of the blade. They depend on 
the equilibritjm bending deflections of the blade, the difference in lead-lag 
and flap and bending stiffness El^ ’ “ Ely? , the torsion frequency and 

the structural coupling parameter (R. For matched stiffness (El 2 » - 
configurations, these couplings vanish. For soft- and stiff- inplane rotor 
blade configurations, they are small and large, respectively. Both torsion 
frequency and precone influence these coupling terms. For very high torsion 
frequencies (o)a > 20) the couplings are normally small and unimportant, but 
for typical values of torsion frequency (5 < < 10) they are important. 

Precone affects the couplings via the equilibrium bending deflections of the 
blades as discussed previously. At low pitch angles, positive precone can 


make Wn 


negative and change the sign of the pitch-lag coupling. For 


0 ' tip 

many rotor blade configurations this will produce instability. Figure 14 shows 



Figure 14.- The effect of torsion frequency and precone on the equivalent 
pitch-lag coupling for zero pitch angle: o)^ = 0.7, (R= 0.0, 0 = 0.0 rad. 
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Figure 17.- Equivalent pitch-flap and 
pitch-lag coupling for a stiff 
inplane rotor blade: = 1.5, 

(R = 0.0, = 0.0 rad. 


Figure 18.- Equivalent pitch-flap and 
pitch-lag coupling for a soft 
inplane rotor blade: = 0.7, 

= 5.0, (R = 0.0. 


General Stability Characteristics, Locus of Roots 

A general overview of hingeless rotor blade flap-lag-torsion stability 
characteristics is next provided by a series of figures showing the locus of 
roots of each fundamental mode as the blade pitch angle 0 is increased from 
zero to 0.5 rad. The main configuration variables of interest are the torsion 
frequency, the lead-lag frequency, and the structural coupling parameter <R. 
Before presenting these results, it will be useful to review the flap-lag 
stability characteristics of a torsionally rigid blade, based on the results of 
reference 6 extended here to include variations in (R (in ref. 6 , (R = 1 only) 
and the second order angle of attack terms. The locus of roots for typical 
soft and stiff inplane configurations are shown in figure 19. 


At zero pitch angle the flap and lead-lag modes are uncoupled. The flap 
mode is highly damped due to the large aerodynamic damping produced by lift 
forces associated with flapping velocity. The lead-lag damping, however, is 
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very small. This is because the aero- 
dynamic drag forces associated with 
lead-lag velocity are very small. In 
the present analysis, structural lead- 
lag damping is not included. Accord- 
ing to results of reference 5, struc- 
tural damping acts like aerodynamic 
profile drag damping, and for weak 
flap-lag instabilities, structural 
damping can eliminate the instability. 

As the blade pitch angle increases 
from zero, the flap and lead-lag modes 
become coupled and the motion may be 
destabilized. The flap mode remains 
stable for both the soft and stiff 
inplane configurations; the lead-lag 
mode is stabilized for the soft 
inplane configuration only, and 
it may become unstable for the stiff 
inplane configuration for values of 
(R between 0 and 1. The effects of 
flap-lag structural coupling are 
larger for the stiff inplane than the 
soft inplane configurations because 
this coupling is proportional to 
El2» - EIy» . 

The flap- lag stability charac- 
teristics will be altered when tor- 
sional flexibility is introduced. 

The changes are expected to be large 
when the torsion frequency is low. 
Figure 20 shows the locus of roots 
for the soft inplane configuration 
for three values of (R and fig- 
ure 21 shows the corresponding locus 
of roots for the stiff inplane config- 
uration. As expected, the effects of 
torsional flexibility are largest 
when C 0 (j) is low and the stiff inplane 
configuration shows larger variations 
in frequency and damping than the 
soft inplane configuration, again due 
to the difference in - EIy». 

The soft inplane lead-lag mode damping 
is generally increased by torsional 
flexibility while the flap mode damp- 
ing is reduced. There is little 
influence of torsion frequency on the 
torsion mode damping. Furthermore, 
the effects of the structural coupling 



(a) Soft inplane, o)^ = 0.7. 



(b) Stiff inplane, = 1.5. 


Figure 19.- Locus of roots for tor- 
sionally rigid rotor blade: 
gpc = 0.0 rad. 
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parameter (R are not large for the soft inplane configuration. The increase 
in lead-lag damping, as torsional flexibility increases, is due to the equiva- 
lent kinematic pitch-lag coupling produced by bending-torsion structural 
coupling. 

The stability of the stiff inplane configuration in figure 21 is strongly 
influenced by torsional flexibility. There is a lead-lag mode instability 
that depends on both and <R and a flap mode divergence for small values of 

(R and Depending on the particular values of (R and o)(f), either the flap 

mode divergence or the lead- lag mode instability will occur at the lowest pitch 
angle. Again, the effects of torsional flexibility can be attributed mainly 
to the bending-torsion structural coupling and the associated kinematic pitch- 
lag coupling. For example, according to the rigid blade flap-lag analysis of 
reference 5, negative 0^ destabilizes the lead-lag mode for small and moder- 
ate values of (R and stabilizes it for (R = 1. For the present results, 9^ 
becomes increasingly negative (fig. 17) as the torsion frequency is decreased 
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•8^9= 0 .3 . 5 , rad 



(c) <R = 1.0. 
Figure 20.- Concluded. 


and, analogous to the rigid blade flap-lag results, the lead-lag mode is 
destabilized when (R = 0 and stabilized when (R = 1.0. At very low torsion 
frequencies, < 2.5, this analog is not valid since torsional dynamics 
terms begin to dominate stability. 


Effect of Lead-Lag Frequency, Torsion Frequency, and Flap-Lag 
Structural Coupling on Stability Boundaries 

In the next series of results, stability boundaries are presented for a 
wide range of parameter values with blade pitch angle 0 the dependent varia- 
ble. At low blade pitch angles, the motion is generally stable, but as 0 is 
increased, the blade motion may become less stable. If this trend continues, 
a stability boundary is crossed at the critical pitch angle 0^^ and the 
motion becomes unstable. In a few cases, the motion is unstable at zero pitch 
angle and may be stabilized at higher pitch angles. Arrows are provided in 
the figures to indicate the stable and unstable regions defined by the 
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stability boundaries and, in some cases, crosshatching is included on the 
unstable side of the boundary. In most cases, the regions of instability are 
clear and the crosshatching is omitted to avoid obscuring the figures. In 
general, an instability of the lead-lag mode is encountered when a stability 
boundary is crossed. In some cases, however, modes other than the lead-lag 
mode determine the stability boundary. In such cases an (F) or (T) , for flap 
mode or torsion mode, respectively, is placed next to the stability boundary. 
The notation ”Div” refers to a flap mode divergence. 

The first result in figure 22 is for the flap-lag stability boundaries of 
a torsionally rigid blade as a function of lead-lag frequency and for several 
values of the structural coupling parameter (R. The soft inplane configura- 
tions are always stable, while most of the stiff inplane configurations 
exhibit a lead- lag instability depending on the values of and (R. These 

pure flap-lag instabilities cannot occur below a certain minimum pitch angle 
which, in this case, is approximately 0.21 rad. As discussed in reference 5, 



(a) (R = 0. 0. 


(b) <R = 0.5. 


Figure 21.- Locus of roots for stiff inplane rotor blade: = 1.5, 

3pc = 0-0 
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the value of this minimum pitch angle depends on the blade profile drag, rotor 
solidity, and blade structural damping (not included in these results). The 
main result to be noted is that as the lead-lag frequency increases, the value 
of (R must also increase to produce instability at a moderate pitch angle. 

For (R > 0.6, flap-lag instability is unlikely for practical lead-lag fre- 
quencies (ci)^ < 2.0). 

Torsional flexibility greatly modifies the flap-lag stability boundaries 
as shown in figures 23 to 28. In each figure, for a single value of (R, 



Wv 


Figure 23.- The effect of torsion frequency on stability boundaries, ver- 
sus lead-lag frequency: = 0,0 rad, (R = 0.0. 



Figure 24.- The effect of torsion frequency on stability boundaries, ver- 
sus lead-lag frequency: Bpc = 0.0 rad, (R = 0.2. 
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Figure 25,- The effect of torsion frequency on the stability boundaries, 0 
versus lead-lag frequency: “ 0.0 rad, (R = 0.4. 



Figure 26.- The effect of torsion frequency on the stability boundaries, 0, 
versus lead-lag frequency: 3r>f> = 0-0 rad, (R = 0.6. 





is varied from 2.5 to In general, soft inplane configurations are free of 

instability for all combinations of (R and o)^. A special limiting case of the 
soft inplane configurations is the matched stiffness configuration where 
Elz' “ y flap-lag and bending-torsion structural coupling vanishes. 

For this configuration the torsion modes are uncoupled from the bending modes, 
and the lead-lag damping may be obtained approximately, but with good accuracy, 
from a flap- lag analysis. Since torsionally rigid soft inplane configurations 
are virtually always stable, torsionally flexible matched stiffness configura- 
tions may be expected to be stable for all values of (R and W(|). This expecta- 
tion is confirmed by the results. Figures 23 and 28 show that the stability 
boundaries shift up and out of the practical pitch angle range as the lead-lag 
frequency decreases to the matched stiffness frequency marked on the figures. 

In contrast to the soft inplane configurations, the stiff inplane config- 
urations again exhibit a variety of unstable behavior. For low values of (R, 
the region of instability occurs at increasingly smaller pitch angles as the 
torsion frequency is reduced. This is directly attributable to the equivalent 
pitch-lag coupling of the bending-torsion structural coupling. As (R 
increases, the region of instability shifts to larger pitch angles and higher 
lead-lag frequencies. When (R = 1.0, the motion is stable for all practical 
configurations. For very low torsion frequency or high lead-lag frequency, 
however, a small unstable region is present. This instability appears to be 
associated with a coincidence of the lead-lag and torsion frequencies 

The results of figures 23-28 are cross-plotted for = 5 in figure 29 
to show more directly the effect of the structural coupling parameter <R. 
Compared to figure 22 for a torsionally rigid blade, the effect of torsional 
flexibility is to reduce significantly the pitch angle at which instability 
occurs. 



Figure 29.- The effect of (R on the stability boundaries, 0^^ versus lead- 
lag frequency: co^ = 5.0, 3p^ = 0.0. 
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In summary, these results without precone indicate that typical hingeless 
rotor blade configurations (o)^ > 5, 0 )^ < 2) operating at moderate pitch angles 
(0 < 0.3 rad) are stable, except when the structural coupling parameter (R is 
relatively small. 


Effect of Precone on Stability Boundaries 

The influence of precone on the stability of torsionally rigid rotor 
blades is relatively small (refs. 5 and 6), but for torsionally flexible 
blades its influence can be very large. This is due primarily to the bending- 
torsion structural coupling terms in the torsion equation. As noted above, 
the equivalent pitch-lag coupling produced by structural coupling is a func- 
tion of the equilibrium flap and lead-lag deflections of the blade that are in 
turn strongly dependent on precone. As clearly shown in figures 14 and 18, 
positive precone produces a positive equivalent pitch-lag coupling, at small 
pitch angles, that can be strongly destabilizing. The stability boundaries 
without precone, given in figures 23-28, are repeated for three values of (R 
in figures 30-32 with 0.05 rad precone. A comparison between the two series 
of figures clearly reveals the destabilizing effects of precone. For low and 
moderate torsion frequencies, a new region of instability appears which begins 
at or near zero pitch angle and terminates at a small positive pitch angle. 

The range of pitch angles for which the blade is unstable increases as 
decreases because the pitch-lag coupling increases with decreasing For 

(R = 0.0 this precone instability is limited to soft inplane configurations; 
however, as (R increases, the region of instability expands to include stiff 
inplane configurations as well. At the low range, near the matched stiff- 

ness configuration, the unstable region closes because the structural coupling 
(EI^^ - ^^y*^ available to generate destabilizing pitch-lag coupling decreases 
and then vanishes when EI^i^ ^^y’ • 
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Figure 30.- The effect of precone and torsion frequency on the stability 
boundaries, 0^^ versus lead-lag frequency: = 0.05 rad, (R = 0.0. 
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In addition to the new region of instability at low pitch angles, precone 
also influences the previous stability boundaries at higher pitch angles, but 
the effect is not very large. Generally, these stability boundaries are 
shifted to higher pitch angles. The effect of increasing precone to 0.1 rad 
for (R = 1 is shown in figure 33; if the results of figures 33 and 32 are 
compared, it will be noted that the regions of instability at low pitch angles 
more than double in size when the precone is doubled. 

Additional results are given for soft and stiff inplane configurations 
where precone is taken to be the independent variable. For the soft inplane 
configuration with o)^ = 0.7 (fig- 34) the motion is stable for > 5.0, and 
only the result for = 2 , 5 is shown. The effect of variations in <R is 



Figure 33.- The effect of precone and torsion frequency on the stability 
boundaries, 0^^- versus lead-lag frequency: =0.1 rad, <R = 1.0. 



Figure 34.- The effect of (R on the stability boundaries, 0^^ versus precone 

angle: 03^ = 2.5, o)^ = 0.7. 
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Therefore, the stability boundary variation with precone for a large torsion 
frequency should be reflected in a similar stability boundary variation over a 
smaller range of precone when the torsion frequency is small. This expecta- 
tion is reasonably well confirmed by the results, particularly for the lower 
values of Also the small effect of precone for large is what would 

be expected from the equivalent pitch-lag coupling effect. 


Effect of Torsion Frequency 

The results discussed above have clearly shown the important influence of 
torsion flexibility. We will now consider this influence more directly by 
showing the variation of stability boundaries as an explicit function of the 
independent variable o)a. These results are shown in figures 39-41 for three 
values of the structural coupling parameter (R and several values of lead-lag 
frequency. For low oi)^, an increase in torsion frequency is stabilizing for 
all values of (R and at higher torsion frequencies, this relationship is 
reversed for some configurations. This occurs, for example, when = 0.0 
near — o)^ and at moderate (R for certain stiff inplane configurations. 

The flap divergence boundaries for low torsion frequencies described earlier 
are shown in detail in the expanded scale portion of figure 40. In general, 
it is difficult to classify the stability characteristics of torsionally flexi- 
ble cantilever blades with respect to the effect of increased or decreased 
torsional flexibility. At large values of the stability boundaries 

asymptotically approach the corresponding flap-lag stability boundaries (pro- 
vided that a flap-lag stability boundary exists) . It may be noted in fig- 
ure 41 that without precone, lead-lag instabilities are not present at typical 
torsion frequencies > 5). Since there are no flap-lag instabilities for 

<R = 1 and 3pe “ 0.0, the results of figure 41 are consistent with the 
results discussed above. 



Figure 39.- The effect of lead- lag frequency on the stability boundaries, 0^^ 
versus torsion frequency: rad, (R = 0.0. 
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Figure 40.- The effect of lead-lag frequency on the stability boundaries, 9^^ 
versus torsion frequency: gp^ = 0.0 rad, (R = 0.4. 

A few results are also included 
to show the effect of torsional flexi- 
bility on lead-lag damping of canti- 
lever rotor blades. This is of inter- 
est because the lead-lag damping of 
cantilever blades is inherently low 
and may lead to problems, such as 
high blade bending loads or coupled 
rotor-body instabilities, even if the 
motion of an isolated blade is stable. 
Figures 42-44 show several examples 
where the damping of stable configura- 
tions is significantly increased or 
decreased by torsion flexibility. 

The damping is given by the dimension- 
less real part of the lead-lag mode 
eigenvalue at a pitch angle of 0.3 rad 
plotted as a function of o)^ for 
? 3 4 ^ several torsion frequencies. The 

influence of torsional flexibility is 
especially clear in these results. 

Figure 41.- The effect of lead-lag As ca^ is reduced, the influence of 

frequency on the stability bound- torsional flexibility is magnified 

aries, 0^.^ versus torsion frequency: whether the lead-lag damping is 

= 0.0 rad, (R = 1.0. increased or decreased. This again 

is due to the influence of bending- 

torsion structural coupling. It is also clear that torsional flexibility 
reverses its Influence when o)^ equals the matched stiffness value. At this 
point, the equivalent pitch-lag coupling changes sign. In general, for 
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typical soft and stiff inplane configurations without precone, torsional 
flexibility increases lead-lag damping except for stiff inplane configurations 
of low (R values. 


Effect of Approximations in the Analysis 

The preceding results and discussion have demonstrated the importance of 
the structural bending- torsion terms in the equations of motion. It is 
natural to inquire, therefore, if the equations of motion could be simplified, 
without significant loss of accuracy, by eliminating most of the terms in the 
torsion equation except the structural terms. This procedure was described in 
a previous section and results in a system of modified flap-lag equations 
without independent torsion degrees of freedom. 

The results using the modified flap- lag equations are compared to the com- 
plete flap-lag-torsion results in figures 45-48. In general, the complete 
flap-lag torsion equations are required for very low torsion frequencies when 
inertial, aerodynamic, and structural terms are all important in the torsion 
equation. At intermediate torsion frequencies, only the structural torsion 
terms are required and the modified flap-lag equations are appropriate. At 
very high torsion frequencies, torsion deflections become negligible and a pure 
flap-lag analysis is adequate. 

An indication of the range of torsion frequency where the three levels of 
complexity are appropriate may be obtained from typical results for a config- 
uration without precone or flap-lag structural coupling ((R = 0) in figure 45. 
For 0 )^ > 20, the flap-lag equations without torsion are adequate; for 
5 < < 20, the modified flap-lag equations including structural torsion 



Figure 45.- The effects of approximations in the analysis on the stability 
boundaries, 0^^ versus torsion frequency: = 0.0 rad, 0.0. 
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terms are sufficient, and for 
( 0 ^ < 5, the full flap-lag-torsion 
equations are required, including 
torsion dynamics. Figure 45 clearly 
shows, however, that this division 
of the torsion frequency range 
depends on other parameters as well. 
For example , with = 1 • 5 , the 
effects of torsion deflections due 
to bending-torsion structural cou- 
pling are significant for > 20. 

Typical results for a configura- 
tion having flap-lag structural cou- 
pling ((R = 1.0, fig. 46) are gener- 
ally consistent with the observations 
based on figure 45. That is, for 
low torsion frequencies < 5, the 
complete equations are necessary, and 
for higher torsion frequencies the 
modified flap-lag equations are 
adequate. 

Figure 47 shows the variation in modal damping versus cda for a stable 
configuration for the three systems of equations. These results are again 
consistent with the range of validity outlined above for the approximate 
equations . 



Figure 47.- The effect of approximations in the analysis on the flap, lead- 
lag, and torsion modal damping versus torsion frequency: = 0.7, 

Bpc = 0.0 rad, 6 == 0.3 rad, <R = 1.0. 



Figure 46.- The effect of approximations 
in the analysis on the stability 
boundaries, 6^,^ versus torsion fre- 
quency: 3pc ” 0-0 tad, (R = 1.0. 
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Figure 48.- The effect of approximations in the analysis on the stability 
boundaries, 9^^ versus torsion frequency: 3p^ “ 0.14 rad, = 1.2, 

(R = 1.0. 

Figure 48 illustrates an example with precone where the modified flap-lag 
equations are accurate to very low for predicting the low pitch angle 

precone instability. A high pitch angle instability at very low is not 

predicted by the approximate equations, however, because of the lack of torsion 
dynamics . 


Effect of the Number and Type of Mode Shapes on Accuracy and Convergence 

The type and number of mode shapes used in the solution of the equations 
has been described previously. It is of interest to examine the influence of 
these factors on the accuracy and especially on the rate of convergence of the 
results. A previous investigation, using uncoupled mode shapes of a nonrotat- 
ing cantilever beam for flap-lag stability analysis (ref. 6), indicated that 
results would be accurate with as few as one or two modes. The present 
results for the flap-lag-torsion equations are not as simple. Figures 49 
and 50 show stability boundaries as a function of precone for a soft and stiff 
inplane configuration, respectively. In each case, the stability boundaries 
are shown for an increasing number of mode shapes for two different calcula- 
tions, the first using uncoupled mode shapes of a nonrotating beam, and the 
second using coupled mode shapes from a free vibration analysis of the rotor 
blade about its equilibrium operating condition (see discussion above in the 
section on Solution of the Equations of Motion) . 

The results using uncoupled mode shapes of a nonrotating beam are given 
in figure 49. These mode shapes yield poor convergence of the stability 
boundaries, and for accurate results five (N) mode shapes are required for 
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0 .02 .04 .06 .08 .10 ,12 .14 


(a) Soft inplane blade: 03^ - 0.7, (b) Stiff inplane blade: co^ = 1.5, 

03^ = 2.5, « = 0.0. o)(() = 5,0, (R = 1.0. 

Figure 49.- Convergence of stability boundaries with respect to number of 
mode shapes for analysis using uncoupled nonrotating mode shapes; total 
number of mode shapes and degrees of freedom equals 3N. 

each bending and torsion deflection. This yields a total of 15 (3N) degrees 
of freedom and is cumbersome computationally. 

The results using coupled mode shapes obtained from a free vibration 
analysis about the equilibrium deflection condition of the rotating beam are 
shown in figure 50. Up to six (M) coupled mode shapes (6 degrees of freedom) 
were used, and for each value of M, 15 (3N) uncoupled nonrotating mode shapes 
were used in the free vibration analysis. It should be noted that a minimum 
of 15 uncoupled nonrotating mode shapes was necessary as a basis for determin- 
ing the coupled rotating modes, this minimum number having been established by 
the results of figure 49. As explained previously, the truncation of the num- 
ber of rotating modes from 15 to M was accomplished in two different ways. 

For M = 5 or 6, the M lowest frequency coupled modes were retained. For 
M = 5 or 6, these modes consisted of 3 or 4 flap bending modes, 1 or 2 lead-lag 
bending modes, and 1 or 2 torsion modes, depending on the particular configura- 
tion parameters. For M = 2, 3, and 4, the modes retained were chosen 



54 



arbitrarily as follows. For M = 4, 2 flap, 1 lead-lag, and 1 torsion mode 
were retained; for M = 3, one of each mode was retained; and for M = 2, the 
lowest frequency flap and lead-lag modes were retained. Figure 50 shows 
clearly that the convergence using coupled rotating mode shapes is quite good; 
in this case as few as three coupled mode shapes provided very accurate results. 


05 \- 


O 


-05 


0 




(a) Soft inplane blade: o)^ = 0.7, (b) Stiff inplane blade: = 1.5, 

03^ = 2.5, (R = 0.0. = 5.0, (R = 1.0. 

Figure 50.- Convergence of stability boundaries with respect to number of mode 
shapes for analysis using coupled rotating mode shapes; total number of 
mode shapes and degrees of freedom equals M. For each value of M, the 
coupled rotating mode shapes are derived from 3N = 15 uncoupled nonrotat- 
ing mode shapes. 


The two widely different configurations used for both figure 49 and fig- 
ure 50 showed similar convergence trends. This indicates that the relative 
effect of using uncoupled nonrotating or coupled rotating mode shapes is not 
highly dependent on the rotor blade configuration. 

In addition to the two ways of varying the number of mode shapes pre- 
sented in figures 49 and 50, it is also possible to use a truncated set of 
coupled rotating mode shapes and vary the number of uncoupled nonrotating mode 
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shapes used as a basis for the rotating inodes. Results were obtained holding 
M constant at 6 with N varied from 2 to 5. Although not shown in the fig- 
ures, the rate of convergence is virtually identical with the results using 
uncoupled rotating modes in figure 49. This is because, for N less than 5, 
the number of nonrotating mode shapes used as a basis for determining the 
rotating mode shapes is insufficient for accurate results. 


CONCLUDING REMARKS 


In this report the stability of elastic cantilever rotor blades of uniform 
mass and stiffness and without twist or chordwise offsets of the mass, elastic, 
and aerodynamic center axes was investigated for the hovering flight condition. 
Instabilities were found mainly for the lead-lag bending degree of freedom. 
These instabilities were strongly dependent on lead-lag frequency, torsion 
frequency, precone, and a structural coupling parameter (R that was used to 
approximate the effects of nonuniformities in bending stiffness. The results 
of the present investigation may be summarized as follows. 

1. The stability characteristics of torsionally flexible blades are 
mainly determined by the nonlinear structural coupling between bending and 
torsion deflections. This coupling produces equivalent kinematic pitch-lag 
and pitch-flap couplings that are generally proportional to the difference 
between lead-lag and flap bending stiffness EI^ » - EI^t , precone, and tor- 
sional flexibility. 

2. Rotor blade configurations without precone, and with typical values 
for torsion and bending frequencies, are generally stable for all practical 
blade pitch angles. Exceptions occur for stiff inplane configurations with 
low values of the structural coupling parameter (R or with low torsional fre- 
quencies, or both. 

3. Rotor blade configurations with precone exhibit lead-lag instability 
for a variety of configuration parameters in practical ranges. This instabil- 
ity exists at low pitch angles and is not present for stiff inplane configura- 
tions when <R is small. 

4. Increasing torsion frequency is usually stabilizing at low values of 

0 )^, but the opposite is true at higher for some configurations. In gen- 

eral, it is difficult to characterize the effects of torsional flexibility in 
a simple manner. 

5. The approximate modified flap-lag equations that neglect torsional 
dynamic effects are usually accurate for torsion frequencies greater than 5.0. 
For 0 )^ < 5.0, torsional dynamics should not be neglected. For torsion fre- 
quencies above 20.0, torsional flexibility is rarely important and flap-lag 
analyses are generally accurate. 

6. The accuracy of the results is dependent on the mode shapes used in 
the modal solution. When using coupled rotating mode shapes, no more than six 
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are required. Using less accurate, uncoupled nonrotating mode shapes, a much 
larger number is required for the same degree of accuracy. 

A logical extension of the present investigation would be to increase the 
scope to include two segment rotor blades with nonuniformities in mass and 
stiffness, twist, and a rigid body blade pitch degree of freedom. Furthermore, 
the effects of droop, sweep, and torque offset, as well as chordwise offsets 
of the elastic, mass, and aerodynamic center axes could be included. 


Ames Research Center 

National Aeronautics and Space Administration 
and 

Ames Directorate, U.S. Army Air Mobility R&D Laboratory 
Moffett Field, California 94035, December 30, 1975 
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APPENDIX A 


DERIVATION OF BLADE VELOCITY COMPONENTS Up, AND h FROM THE 

DEFORMED BLADE COORDINATE TRANSFORMATION 


In equations (18) and (20) blade aerod 3 mamic loads are expressed in terms 
of U-j;, Up, and e where Up and Up are components of blade velocity along 
the y’ and z* axes, and e is the component of blade angular velocity along 
the x’ axis. In order to use the expressions for blade aerodynamic loads in 
equations (10), Up, Up, and e must be expressed in terms of the blade bending 
and torsion deflections v, w, and <p. The blade velocity is easily expressed 
in the x, y, z coordinate system. The deformed blade coordinate transforma- 
tion [T] relating the x, y, z and x’ , y’, z’ coordinate systems, derived 
in references 1 and 24, is then applied to express the Up and Up velocity 
components in terms of v, w, and (p. The angular velocity components of the 
blade may be obtained by taking advantage of certain mathematical properties 
of [T]. 

The vector velocity U of the blade may be simply written in the x, y, 
z coordinate system from equation (45) of reference 1, and with the addition 
of uniform induced inflow vp becomes 

U = (u - ^^v)i + (Qx + v) j + (vp + w + f^vBp^)^ = ^R^' ^T^ * Upk* 


The transformation [T] relates the k, and i’ , j ’ , k’ unit vectors 

of the X, y, z and x’ , y’ , z’ coordinate systems, respectively. There- 
fore, the blade velocity components in the deformed x’ , y’ , z’ coordinate 
system are 


Ur 


' u - Qv 

Uj 

= [T] 

Qx + V 

Up] 


i + w + j 


(Al) 


where, valid to second order. 


[T] = 


1 - 




•[v' cos(0 + (|)) + w’ sin(6 + c^))] 


[v* sin(6 + 4)) - cos(0 + 4>)] 


-(■ - 


cos (0 + } + ) 

in (0 + $ + ) 


(' ■ " 0 " 


sin(6 + <^) 
cos(0 + 4 >) 


and 


$ = (p 


f v"w' dx 
0 


(A2) 
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Thus, using equations (Al) and (A2) 

U.JT = - (u - ^2v) [v’ cos (0 + cj)) + w’ sin (0 + (j>) + (to + v) cos (0 + ^ + v’w') 

+ (vi + w + to3p^)sin(0 + $) 

Up = (u - to) [v’ sin(0 + (j)) - w' cos (0 + (f)) ] - (to + v) sin(0 + $ + v’w’)^l - 
+ (vj^ + w + to3p^)cos(0 + $)^1 - 

For 0 = 0(e) we have sin 0 = 0 and cos 0=1. With all terms of O(e^) 
neglected with respect to unity, the and Up velocity components are 

U-j = Uk + V 

X \ (A3) 

v’w" dx I + v^ + w - (0 + (p)v + to(3p(, + w’) 




The component of angular velocity along the x' axis, e, is now con- 
sidered. First, an infinitesimal rotation [a] dt is defined where 


[a] 


0 a 


“Ot. 


^ J 

0 a . 


a . -a . 0 

3 1 


(A4) 


and where aj , and are the angular velocity components about i’ , j’, 

and k’ . Note that e = a;j^. By definition, the matrix [T] satisfies the 
relation 



Differentiating this equation yields 


dt- 

1^1 

dj' 

II 

H 


-> 

dk’ , 

Ikl 


which in turn can be written in terms of [a] as 
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Thus, 


[a]dt<j'V = [dT]<t 


[T] = [a][T] 


(A5) 


The transformation [T] was originally derived in reference 24 by solving the 
following differential equation, similar to equation (A5) : 

[T'l = M[T] (A6) 

where [cj] is of the same form as [a] and the elements wi, wj , to^ are 
rates of rotation with respect to the blade axial distance x, functions of 
v", w", and (J)'. We now consider the derivatives of the tensor function 
[T(v' ,w' ,(())]: 


[T-] = 
[T] = 


8T „ . 8T „ ^ 9T 

3v' 3w' ^ 9(|) 


3T . , . 3T . , . 9T : 


(A7) 


By comparing the two equations of (A7) and equations (A5) and (A6) it is 
evident that, when v" , w” , and (f) ’ in the expressions for 03 j , and co]^ are 

replaced by v’ , w’ , and the exact expressions for a^, a j , and a]^ are 
obtained for a nonrotating blade. Making this replacement in the expressions 
for 0 ) j , and in the appendix of reference 1 and adding the i’ , j’, k' 

components of the blade angular velocity ^ of a rotating blade yields 



v’ sin(0 + cj)) - w 
v' cos(0 + (())+ w' 


cos(0 + (\>) 
sin(0 + (j)) 


+ a[T] 



(A8) 


The only component of interest here is = e. Thus, neglecting higher order 

terms 


e = + w’) 


(A9) 
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APPENDIX B 


MODAL EQUATIONS OF MOTION 


In this appendix the nonlinear, variable coefficient, integro-partial 
differential equations (10) and (23) are written in nondimensional form and 
transformed into nonlinear ordinary differential equations. These modal equa- 
tions are then linearized for small perturbation motions about the equilibrium 
operating condition. Nonlinear algebraic equations, governing the equilibrium 
operating condition, are written first, and the mass, damping, and stiffness 
matrices, equation (28), for the linearized perturbation equation are then 
expressed in terms of the equilibrium generalized coordinates. 

First, equations (10) and equations (23) for the aerodynamic loadings are 
combined and written in nondimensional form: 


6v equation: 

.1 


+ A 2 V"" - (A2 - Apsin2((Re)v"" + - A^ w"" 


v' J* (x + 2v)dx 

X 

+ (A 2 - Ap [-sin(2(R0) ((|)v")” + cos (2<R0) ((J)w”)” ] - 23p^.w - 2 J (v’ v’ + w’ w’ ) dx 


+ V - V + ^ 


XV£(f) + 


2 X + (0 + (|))v. 

a 1 


J 


V + [2v^ - x(0 + (j))]w 


6w equation: 
.1 


= ^ ^ x2 - xviej 

(Bla) 


w’ / (x + 2?)dx ' + Ajw"" + (A 2 - A^)sin2((Re)w"" + (A 2 ~ A^) v’’" 

X J 

+ (A 2 - Aj^) [cos(2(Re) ((j)v")" + sin(2(Re) ((J)w")"] + 23p^v + w ^1 + ^ j 

-x2^(j) + J v'w" djc'j + xv(3pQ + w* ) - xw' - [2x(6 + <|)) - v^]v 




3c — • 
+ xw - X (|) 


= -^pc^ + 6 


+ X^e + xBpj,^ (Bib) 
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6<1) equation: 




[(1 - x2)(j)']' - Kf + 
— 2 


<h - A,)[(5!^if^)sln(2«6) 


+ v”w" cos(2(R0) 


'] 


+ x4> + + (\i2^ ~ cos(20) = “(y 2 ^ “ Vii^) (Blc) 
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The dimensionless parameters are given by 


K = 


m 


‘m. 


- X 
X =- 


V = 


V 

R 


Ml = 


1 R 


An = 


Ely, 


ml^^R 


2p4 


_ w 
^ = R 


■ 


2 R 


R 


_ c 
^ = R 


_ EI^’ 




^1 J2R 


GJ 

^|!S 

11 

) 


3p acR 

(•) = — = 
'■ ^ 9i|) 

3 

m 

J23t 


( ) 


(B2) 


Substitution of equation (25) into equations (Bl) and application of Galerkin’s 
method yields the following set of nonlinear ordinary differential equations: 


N 


6V^ equations: ^ | Vj{D^j + [A^ - (A^ - A^ ) sin^ ((R0) ] g j ^6 - ^ij ^ 


J = 1 


+ Wj (A^ - + (A 2 - A^) 2 Kjki^j[-Vk sln(2(Re) 

k — X 

N ^ N , 

+ cos(2(R0)] - 23pc5ijWj + 2 (^ikj “ ^ 

k— 1 X— 1 

/ N \ 

'I 0 + (2 ^ E.. + 0V.6.. + V. y )v. 

1 Ij j \ a ij X ij ^ ^ 


+ 6 . .V. + -z 
ij J 6 


I ^ \- 

+ 2v.6.. - 0E.. - y $-G..,)w 

\ 1 ij iJ 1^1 ^ 


= ^ ( V . ^A. C . - V . 6B . j 1-1,2, . . .,N 

V ^ ^ ^ ^ V (B3a) 
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6W^ equations: 

^ ' sin(2(Re) 




N 


+ (A 2 - Ai)sin2((Re)]3 /"a } + (A^ - A^ ^ ^ikl^ cos(2fl0) 

J J ^ 

+ .in(2fie)] + 2Bp,JyVj + 2 £ + ii^ (l + 




N 

“j..^. + ]T L.mV.W. + S„pE.,V. - f 0..W. 
IJ 3 ijl<: 3 J IJ J 2 XJ j 


N 


- 20E 


3c 


+ 2 52 4’iG,., -v,6..)V. + E..W. - -I..0. 
13 k xjk X xj/ j ij j 4 xj 3 


■ -®P='i + 6 (-Vl + “<=1 + f »pc®i) 


1=1,2 N 


6$^ equations: 
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2 ( $ (m^kn^ + [ky 2 + - y^2)£;o3(2e)](S. . } 

j — \ J J J -■- J 


+ (A 2 - A^ ^ K 
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ijk 


+ cos(2(Re)V^W^ 


+ . .$. + M. .$. 

ij J 48 xj J, 


= -(y^^ - y^2) slr^ 


/2Y^ 


1=1,2, 


(B3b) 


N (B3c) 
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Here 


L = f 

X < 1 


. . = f ^ x^'i'.e. 

ij •'o 1 1 


" dx 
k 


:. = / x^f. 

X < X 


1 

L. ., = f x'i'.f.'i', ' dx 

xjk J X j k 


-77^ h' dx 

2 / X j 
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13 Jo 1 
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The equilibrium and perturbation equations are obtained by substituting 
equations (27) into equations (B3) and performing the operations described in 
the text. The resulting equilibrium equations are: 


^ equations: 2 [^2 “ (^2 ~ Ai)sin2((R6)]3j^6Fj - 


sin(2(R6) 


Oj 2 


l2-Al)6.-«y + (A^-Ai) E sin(2fie) 

k 1 


+ cos(2<R0)] + ^ Vi/oj)= i I 


— C. - V.0B. i= 1, 2, . . . , N 

^ 7 (B5a) 
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6W^ . equations: 
Ox 


E (v„, (A, - . [A, 

N 

+ (A 2 - Apsln2((Re)]B^^6^^} + (A^ - A^ ^ ^jkl^Oj ^°s(2(R0) 

+ sin(2(Re)] 
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6 \"Pc-ij^0j 2 -ij^^Oj -ij-Oj 
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6 $ 0 ^ equat ions : 


2 ■•■ (^2^ " y;^^)cos(26)]6^_ } 
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+ (A^-Aj) 2 K 


k=l 


ijk 


iJ 


«Oj"0k - 'o/ok) 1- “*<2««''oj''ok 


- -(P,2 - P.2) 1-1,2,. 

/2Y^ 


, N (B5b) 


, N (B5c) 


The perturbation equations, linearized in A-quantities , yield a set of 
ordinary differential equations with constant coefficients that can be easily 
expressed in the form of equation (28) where 



(B6) 
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